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ABSTRACT 

This paper illustrates how mock observational samples of high-redshift galaxies with 
sophisticated selection criteria can be extracted from the predictions of GALICS, a 
hybrid model of hierarchical galaxy formation that couples the outputs of large cos- 
1"***" . mological simulations, and semi-analytic recipes to describe dark matter collapse, and 

' the physics of baryons respectively. As an example of this method, we focus on the 

properties of Lyman Break Galaxies at redshift z ~ 3 (hereafter LBGs) in a ACDM 
cosmology. With the MOMAF software package described in a companion paper, we 
generate a mock observational sample with selection criteria as similar as possib le to 
those implied in the actual observations of z ~ 3 LBGs bv ISteidel et al.l l)1995|h We 
need to introduce an additional "maturity" criterion to circonvene subtle effects due 
to mass resolution in the simulation. We predict a number density of 1.15 arcmin 2 at 
q R < 25.5, in good agreement with the observed number density 1.2 ±0.18 arcmin -2 . 

5_i ■ Our model allows us to study the efficiency of the selection criterion to capture z ~ 3 

^ \ galaxies. We find that the colour contours designed from models of spectrophotomet- 

er ■ ric evolution of stellar populations are able to select more "realistic" galaxies issued 

from models of hierarchical galaxy formation. We quantify the fraction of interlopers 
(12 %), and the selection efficiency (85%), and we give estimates of the cosmic vari- 
ance. We then study the clustering properties of our model LBGs. They are hosted 
by halos with masses ~ 1.6 x 1O 12 M0, with a linear bias parameter that decreases 
with increassing scale from b — 5 to 3. The amplitude and slope of the 2D correlation 
function is in good agreement with the data. We investigate a series of physical prop- 
erties: UV extinction (a typical factor 6.2 at 1600 A), stellar masses, metallicities, and 
Star Formation Rates, and we find them to be in general agreement with observed 
values. The model also allows us to make predictions at other optical and IR/submm 
wavelengths, that are easily accessible though queries to a web interfaced relational 
database. Looking into the future of these LBGs, we predict that 75 % of them end 
up as massive ellipticals and lenticulars today, even though only 35 % of all our local 
ellipticals and lenticulars are predicted to have a LBG progenitor. In spite of some 
shortcomings that come from our simplifying assumptions and the subtle propagation 
of mass resolution effects, this new 'mock observation' method clearly represents a 
first step toward a more accurate comparison between hierarchical models of galaxy 
formation and real observational surveys. 

Key words: galaxies:high-redshift - galaxies:formation - galaxies:evolution - astro- 
nomical data bases:miscellaneous 
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1 INTRODUCTION 

Models of hierarchical galaxy formation are now sophisti- 
cated enough to produce a host of predictions for the statis- 
tical properties of local and high-redshift galaxies. These 
models have to be tested against observational samples. 



However, such a comparison is not as straightforward as 
it might appear at first glance, because observations are af- 
fected by selection criteria and observational biases. The dif- 
ficulty of the comparison is enhanced for samples of high- 
redshift galaxies that are selected only on the basis of their 
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apparent magnitudes and colours, that is, according to prop- 
erties that are far from the dynamical quantities models of 
galaxy formation deal with. In such a case, it is not easy to 
track back the propagation of the selection criteria and ob- 
servational biases to the "model" parameter space. At this 
stage, the best method consists in putting model predictions 
into the "observation" parameter space, by producing mock 
galaxy samples that are obtained with selection criteria and 
observational biases as close as possible to those that affect 
the real observational samples. 

This paper uses the predictions of the GALICS model of 
hierarchical galaxy formation (for Galaxies in Cosmological 
Simulations) that embodies the so-called "hybrid approach" 
by coupling the description of dark matter collapse by means 
of cosmological iV-body simulations, and the description of 
the physics of baryons with semi-analytic recipes, including 
the estimate of spectral energy distributions in the UV to 
sub- millimetre wave length range jHafrton et al.l(l2003T) .here- 
after GALICS I, and iDevriendt et aTn200^l . hereafter GAL- 
ICS n). We use the momaf packa ge (for Mock Map Facility) 
described in lBlaizot et alJ <|2003T) (hereafter momaf) to pro- 
duce mock observing cones from which mock samples can 
be extracted, with selection criteria and biases that mimic 
those of actual observations. 

We hereafter illustrate such a method by addressing the 
constraints put on models by the so-called Lyman Break 
Gala xies (hereafter LBGs) at redshift z ~ 3 iSteidel et alJ 
1996). These galaxies are obtained using the UV drop- 
out techniq ue, first employed f o r the selection of distant 
galaxies by ISteidel &z Hamilton! (|l993j), that simply relies 
on the shift of the Lyman break through broad-band fil- 
ters. This break, being mainly due to the absorption of 
the UV photons by the Hi of the intergalactic medium 
(IGM), is roughly independent from the galaxies' intrinsic 
properties, thus allowing the efficient selection of a com- 
plete sample of luminous galaxies at a given redshift, with 
photometry in only three broad-band filters. In the last 
years, this technique has led to a tremendous increase in 
the number of observed objects at z ~ 3. Although the 
set of data is impressively large today, compared to only 
a few years ago, a theoretical picture has not yet clearly 
emerged. This is mainly due to a poor understanding of 
the interplay of star formation, feedback and extinction 
on cosmological scales within hierarchical clustering, and 
also to the difficulty to compare "unbiased" and "noiseless" 
theoretical work with "corrected-for-everything" data. In 
spite of these difficulties, the predictions of the properties 
of LBGs obtained fro m semi-analytic m o dels of hierarchi- 
cal ga la xy formati on fe augh et all <ll998l) , ISomerville et al] 
{200 ll). Edzi et al.N2003h1 or fully num erical approaches 
jNagamind foOOfl . I Weinberg et all j2002D l are encouraging. 

We hereafter produce a mock sample of z ~ 3 LBGs 
from our model of hierarchical galaxy formation, by apply- 
ing selection criteria and observational biases that are as 
similar as possible to those invo l ved in the definition of the 
data gathered by ISteidel et alJ il99rJ) . The comparison of 
our mock sample with the actual sample has a three-fold 
interest. First, it helps us to explore the selection criteria 
with a plausible model of galaxy formation, whereas the UV- 
drop out technique has been primarily designed within the 
paradigm of monolithic galaxy formation. This allows us to 
get some insight on the efficiency of the selection criteria to 



capture bona fide galaxies at z ~ 3, as well as to assess the 
number of interlopers, that is, of galaxies that pass through 
the criteria, but are not at the target redshift z ~ 3. In addi- 
tion we study the field-to-field cosmic variance. Second, we 
compare predictions of our model of hierarchical galaxy for- 
mation with observational properties. Since GALICS has spe- 
cific features such as spatial information, multi-wavelength 
predictions from the (rest-frame) UV to sub-mm, and hier- 
archical information on the merging history trees, we can 
address questions such as clustering, optical/IR luminosity 
budget, and current descendants of z ~ 3 LBGs. Third, our 
predictions can be used as an educated guideline for exten- 
sive follow-up observations of LBGs at other optical and 
IR/submm wavelengths. The overall agreement of the pre- 
dictions with the data, in spite of obvious shortcomings, en- 
courages us to think that the picture of hierarchical galaxy 
formation implemented in our GALICS model gives a good de- 
scription of optically-bright, star forming galaxies at z ~ 3. 
We note however that not all galaxies at this redshift are 
correctly reproduced by GALICS. 

This paper is the third in the GALICS series after pa- 
per I (GALICS i) that presented the global features of our 
model and a set of properties for local galaxies, and paper 
II (GALICS ii) that showed predictions for hierarchical galaxy 
evolution between redshifts and 3. It is also an illustration 
of the generic method of mock map making detailed in mo- 
maf. Predictions of faint galaxy counts and 2D clustering 
for the global population will be given in forthcoming paper 

IV (Devriendt et al. 2003b, hereafter GALICS iv), and paper 

V (Blaizot et al. 2003b, hereafter GALICS v) respectively. 
These five GALICS papers, together with the momaf pa- 
per, complete the presentation of the first version of the 
GALICS project. Results of the model, as well as extensive 
follow-up predictions at many optical and IR/submm wave- 
lengths, are publicly available from a web interfaced rela- 
tional database located at the URL http : //galics . iap . f r 

This paper is organised as follows. Section 2 gives a brief 
summary of the features of GALICS and MOMAF that are rel- 
evant for this study. The selection criteria are applied to our 
mock observing cone in section 3, where we explore various 
aspects of the UV-drop out technique. Section 4 gives pre- 
dictions for the cosmic variance, and compares the predicted 
clustering with observational data. The UV to IR luminos- 
ity budget is studied in section 5, where we try to clarify 
the controversial issue of extinction in LBGs from a theo- 
retical point of view. An attempt to elucidate the nature of 
LBGs, and their fate at the present time within the context 
of hierarchical clustering, is given in section 6. Finally, sec- 
tion 7 gives a summary of the results, and discusses their 
robustness in light of the shortcomings of the method. 



2 SIMULATING OBSERVATIONS 

2.1 A brief summary of the GALICS model 

GALICS is a model of hierarchical galaxy formation which 
combines high-resolution cosmological simulations to de- 
scribe the dark matter content of the Universe with semi- 
analytic prescriptions to deal with the baryonic matter. This 
hybrid approach is fully described in GALICS I and GALICS II, 
so we only briefly recall its main characteristics hereafter. 
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2.1.1 Dark matter simulation 

The cosmological N-body simulation we refer to through- 
out this paper w as carr i ed ou t using the parallel tree- 
code developed bv iNininl jl999T) . We use a flat Cold Dark 
Matter model with a cosmological constant (Q m = 0.333, 
f^A = 0.667). The simulated volume is a cube of side 
Lbox = 100/i^QoMpc, with h^ — 0.667, containing 256 3 
particles of mass 8.3 x 1O 9 M0, with a smoothing length of 
29.29 kpc. The power spectrum was set in agreement with 
the present day abundance of rich clusters (as = 0.88), and 
we followed the DM density field from z = 35.59 to z = 0, 
outputting 100 snapshots spaced logarithmically in the ex- 
pansion factor. 

On each snapshot we use a friend-of-friend algorithm 
to identify virialised groups of more than 20 particles, thus 
setting the minimum dark matter halo mass to 1.66 x 10 
Mq. We compute a set of properties of these haloes, includ- 
ing position and velocity of the centre of mass, kinetic and 
potential energies, and spin parameter. Assuming a singu- 
lar isothermal sphere density profile for our virialised dark 
matter haloes, we then compute their virial radius, enforcing 
conservation of mass and energy. 

Once all the haloes are identified in each snapshot, 
we compute their merging history trees, following the con- 
stituent particles from one output to the next. The merging 
histories we obtain are by far more complex than in semi- 
analytic approaches as they include evaporation and frag- 
mentation of haloes. The way we deal with these is described 
in details in GALICS I. 



2.1.2 Baryonic prescriptions 

When a halo is first identified, it is assigned a mass of hot 
gas, assuming a universal baryon to dark matter mass ratio 
(fis/fio = 0.135 in our fiducial model). This hot gas is as- 
sumed to be shock heated to the virial temperature of the 
halo, and in hydrostatic equilibrium within the dark matter 
potential well. The comparison of the cooling time of this 
gas to its free-fall time, as a function of the radius, yields 
the mass of gas that can cool to a central disc during a 
time-step. The size of this exponential disc is given by con- 
servation of specific angular momentum during the gas infall 
and scales linearly with the spin parameter of the halo. Cold 
gas is then transformed into stars at a rate inversely propor- 
tional to the dynamical timescale, with a constant efficiency 
0~ . Newly form ed stars are distributed according to the 
iKennicuttl il983t) initial mass function (IMF). These stars 
are then evolved between timesteps, using a sub-stepping 
of at most 1 Myr. During each sub-step, a fraction of the 
stars explode as supernovae, releasing metals into the ISM. 
Energy is also released, which heats the ISM and may blow 
part of it away into the IGM, with some efficiency e. 

When two haloes merge, the galaxies they contain are 
placed in the descendant halo. Their initial radial positions 
in this new halo are obtained through a perturbation of 
the final radial positions they had in each of the progenitor 
haloes. Due to subsequent dynamical friction or satellite- 
satellite encounters, galaxies can then merge in their new 
host halo. When such events occur, a 'new' galaxy is then 
created (which is the descendant of the two progenitor galax- 
ies that have merged) and the stellar and gaseous contents 



of its three components (disc, bulge and starburst) are de- 
duced from those of its two progenitors using a single free 
parameter \. Note that the descendant galaxy can either 
become elliptical (in shape) or retain a large disc depending 
on the mass ratio of its progenitors. After a merging, pro- 
vided gas is available, a nuclear starburst is assumed to take 
place at the centre of the 'new' galaxy. The radius of this 
burst is fixed to be k times the radius of the 'new' bulge. 

The spectral energy distributions (SEDs) of our mod- 
elled galaxies are computed by summing the contributions of 
all the stars they contain, according to their age and metal- 
licity, both of which we keep track of in the simulation. Ex- 
tinction is computed assuming a random inclination for disc 
components, and the emission of dust is added to the ex- 
tinguished stellar sp ectra in a similar way as prescribed in 
the STARDUST code iDevriendt et aljll999f) . Finally, a mean 
intergalactic medium (IGM) extinction corr ection is imple- 
mented following the procedure described in iMadaul lll995l) . 
before we convolve the SEDs with the desired observer frame 
filters (see momaf). 

Thus our model depends only on a few free parameters, 
which are given the standard values quoted in GALICS I and 
GALICS II, namely, /3 _1 = 0.02, e = 0.1, % = 3-333, n = 0.1. 
As it is explained in GALICS I, the other two free parameters 
(the normalisation ip of the frequency of satellite-satellite 
encounters, and the efficiency £ for the recycling of metals 
ejected into the intergalactic "reservoir" ) only have a minor 
role. 



2.2 Resolution effects 

The mass resolution of the N-body DM simulation we use 
has important repercussions on the statistical and physical 
properties of our modelled galaxies. We outline the main two 
effects below. 



2.2.1 Magnitude completeness limit 

The mass resolution of the N-body simulation is, for our 
concerns, the minimum mass of a halo, that is, the mass of 
20 particles (~ 1.6 x 10 U M Q ). It is the mass of the smallest 
structure in which we may form a galaxy. This halo mass 
corresponds to a galaxy mass M res = 2.2xlO lo M assuming 
that all the gas in the halo can cool. Galaxies with a mass 
higher than this cannot exist in unresolved haloes, and our 
sample of such galaxies is thus complete. On the contrary, 
galaxies with a lower mass can (and do) exist, but our sample 
is not complete since we miss all those that would lie in 
unresolved haloes. 

It is useful to convert this completeness limit in terms 
of magnitudes in order to understand in which range of lu- 
minosity our predictions are reliable (e.g. for the luminosity 
functions in section 15. H . To do this, we plot the absolute 
rest-frame UV magnitude at 1600A (internal dust absorp- 
tion included) versus the total baryonic mass of galaxies in 
Fig. Q directly from the z — 3 simulation snapshot. Note 
that the large scatter on this plot is due to the fact that at 
this wavelength, we see the combined effects of most recent 
star formation and dust extinction, which are poorly corre- 
lated with the total mass of stars. We define the magnitude 
completeness limit M^£(1600) so that 90% of the galaxies 
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Figure 1. Absolute AB magnitudes of z = 3 galaxies at rest- 
frame 1600 A as a function of baryonic mass. The absolute mag- 
nitudes take into account dust absorption. The vertical line indi- 



cates our mass resolution, M r 



2.2 1O 1U M , above which our 



sample of galaxies is complete. The horizontal line represents our 
completeness limit in terms of magnitude so that galaxies brighter 
than M^"^(1600) = —20.3 also form a complete sample. 



brighter than M^(1600) have a mass greater than M res . 
As shown in Fig.Q we expect our sample of galaxies brighter 
than Mis (1600) = -20.3 to be complete at z ~ 3. 

The magnitude completeness limit of our simulation is 
thus comparable to the detection limit for Steidel's sample 
of LBGs (absolute magnitude R < -20.2 in our ACDM 
cosmology), which allows us to compare our results with 
their data. 



2.2.2 History resolution 

A more subtle effect of resolution is that missing small struc- 
tures means missing part of our halo merging histories. As 
a matter of fact, the impact on our modelled galaxies is 
twofold : 

• we miss small galaxy mergers because galaxies close to 
the resolution limit are formed smoothly in haloes at the 
resolution limit, instead of assembled through mergers of 
smaller objects that would have populated smaller progeni- 
tor haloes; 

• we cool too much gas in too short a time on new central 
galaxies in haloes at the resolution limit. Had we resolved 
their host halo merging histories, part of this cold gas would 
have been split between smaller progenitors at earlier times 
and partly processed into stars already. 

However, the 'observational' signature of this resolution ef- 
fect (physical properties, luminosities of galaxies) tends to 
disappear after a while. In practice, a galaxy in our standard 
N-body simulation needs to evolve for about 1 Gyr before 
we can be sure that its properties have converged. This is 
irrelevant at z = (resp. z — 1), where only 3 (resp. 69) out 
of 19,372 (resp. 19,257) galaxies above the formal resolution 
limit Ms = -18.9 (resp. Ms = -20.2) are younger than 
this threshold. However, it becomes problematic at z = 3, 
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Figure 2. Mock medium deep field, of size 3x3 arcmin 2 , seen 
in the R band. Circles mark the position of LBGs which were 
selected using the colour and apparent magnitude criteria given 
by Eq.0 



when the age of the Universe is only a couple of gigayears 
old and roughly 75% of the galaxies are concerned. As a con- 
sequence, we introduce in the following study an additional 
free parameter that we call the maturity age t ma t, and we 
discard from our sample any LBG whose first progenitor is 
younger than t m at = 1.1 Gyr at the time when this LBG is 
identified as such. As the need for this kind of selection is one 
of the most important drawbacks of our study, section 3.2.2 
justifies in detail our choice of maturity criterion. 



2.3 A typical mock field 

Because Lyman Break selection is by definition an observa- 
tional selection, it is necessary to apply as similar as possible 
criteria to our modelled galaxies sample before attempting a 
thorough comparison with the data. We therefore used the 
MOMAF package to simulate light-cones from GALICS out- 
puts, as described in the MOMAF paper (see also the GAL- 
ics/momaf web-page http://galics.iap.fr/]. Basically, 
the method consists in tiling a time sequence of simulation 
boxes along a mock line of sight, and computing the appar- 
ent properties of galaxies from their positions in this cone- 
like geometry. It allows us to use as much of the spatial 
information contained in the N-body simulation as possi- 
ble in order to produce realistic mock catalogues, which can 
then be analysed in the same way as real observations. 

In Fig. |5J we show a typical result of this work : a 
mock medium deep field of 3 x 3 arcmin 2 seen in the R 
band ijSteidel fc Hamiltonll993F) . The circles indicate LBGs, 
which are identified as indicated in section |3~T1 
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3 LBG SELECTION PROCESS 
3.1 Colour-colour selection 

It is the IGM extinction that makes the Lyman break se- 
lection such an efficient technique for selecting high-redshift 
galaxies. Galaxy evolution, along with dust absorption in- 
side galaxies, also produce a Lyman-break, but its intensity 
varies from one galaxy to the next, so that galaxies end up 
scattered all over the colour-colour plane and as a result, 
colour and redshift would be strongly degenerate. Not only 
does the IGM attenuation yield a more pronounced Lyman- 
break, but the fact that it is due to material external to the 
galaxies makes this break independent from intrinsic prop- 
erties of the observed objects, allowing the observer to select 
a whole population of "normal" galaxies at z ~ 3. 

For a galaxy at this redsh ift, the Lyman break fall s 
between the U„ and G filters (Stejdel fc Ha milton! Il993l) . 
It makes the U„ — G colour of galaxies very sensitive to 
small variations of redshift around this value, as shown by 
the very steep slope in the colour-colour diagram of Fig. |3] 
This spread of the galaxy distribution along the U„ — G axis 
leads to a precise selection of galaxies in redshift space. The 
colour criteria for selecting objects at z ~ 3 proposed by 
ISteidel et alJ l|l995|l . which we adopt throughout this paper, 
is the following : 
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In Fig. |3] we show a colour-colour diagram for a mock 
catalogue of 360 arcmin 2 . Only galaxies with an apparent 
AB magnitude lower than 25.5 in the R band are shown, in 
order to mimic the detection limit quoted by these authors. 
We represent galaxies which have redshifts between 2.7 and 
3.4 as diamonds. Galaxies with other redshifts are marked 
as simple dots. The colour criterion given by Eq.^is shown 
as the solid line. 

Notice that the distribution of galaxies in the colour- 
colour plane is smooth in Fig. |3] This behaviour stems 
from distributing galaxies in a mock light-cone, in contrast 
with previous studi e s based on pure s e mi-an alytical models 
jBaueh et alJ (ll998Tl . ISomerville et alJ (1200 lM . where galaxy 
properties were only computed for a discrete set of time out- 
puts. 

3.2 Data sets 

3.2.1 Observations 

We choose to compare our modelled LBGs to the sample 
presented by Steidel and co-workers, for both the variety 
of res u lts they have publis h ed so far (e.g. iGiavalisco et alJ 
<199Sl) . ISteidel et al] <1999MShaplev et alJ <200lft ). and the 
large number of sources it contains. The total area surveyed 
by these authors amounts to about a th ousand arcmin 2 , 
and co ntains ~ 1250 LBG candidates (see iGiavalisco et alJ 
( 1998), table 1, for a summary of observations). The density 
of observed LBGs on the sky is thus about 1.2 object per 
square arcminute, among which roughly 20% are estimated 
to be interlopers. Most of the observations were made using 
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Figure 3. (Un — G)-(G — R) diagram, for a mock 360 arcmin 2 
field. Only galaxies with apparent magnitude Rab < 25.5 are 
shown. Diamonds represent galaxies that have their redshift in 
the range [2.7; 3.4], and points show galaxies out of this range. 
The solid contour shows the colour criteria for LBGs at z ~ 3 
(Eqs. , and the dashed contour shows our optimised contour 
for selecting galaxies at z ~ 3 fsection !3.4l . Notice the very steep 
rise in U n — G around z ~ 3, mainly due to the IGM attenuation. 
Photometric errors are not added in this plot. 



the three filters U n , G, and R, designed to efficiently select 
objects at z ~ 3 ISteidel fc Hamiltonlll993T) . 

A point which is still unclear in the literature so far is 
the definition of interlopers. In principle, these are galaxies 
which, although selected by the colour criterion, are not at 
the requested redshift. In practise, it seems that the defini- 
tion which is used by some authors is somewhat less strict, 
and objects are only called interlopers when they are galax- 
ies at redshifts lower than ~ 2.3 or stars. We follow suit and 
adopt this broader definition but remark that since there 
are no foreground stars in our simulation, we do not expect 
to find many interlopers (see following section). 

3.2.2 Simulations 

We us e the (U n , G, R) filter set from ISteidel fc Hamilton! 
(1993) to compute the apparent magnitudes of our simu- 
lated galaxies. Our standard mock catalogue size is 1 deg 2 , 
and contains sources up to z ~ 6.5. We also use the same 
apparent magnitude detection limit, that is, R < 25.5. 

As we discuss in momaf, a unique way to build a mock 
catalogue from the outputs of our simulation does not ex- 
ist, because the mock light-cone intersects only fractions of 
each snapshot. We choose to introduce a random process in 
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Figu re 4. LBG count s from GALICS, compared to data from 
iGiavalisco &: Dickinson! 1200 J) (diamonds with error bars). The 
solid line and grey area show the mean and l-a region obtained 
from 20 mock catalogues of 1 square degree in which we applied 
the maturity age selection with t ma t = l.lGyr. The dashed line 
shows the LBG number counts for a single cone, without any ma- 
turity selection, and the dot-dashed line shows the same counts 
with t ma t = 1.3 Gyr. The dotted line shows counts from a higher 
resolution simulation covering a much smaller volume. Its normal- 
isation justifies the value given to tmat (see text for discussion). 

the generation of cones in order to be able to make inde- 
pendent realisations which are all a priori equally plausible. 
Because we do not include fluctuations on scales larger than 
the simulated volume, which is about 1 degree at z ~ 3, the 
dispersion of the number of LBGs found in different cones 
will yield an under-estimate of the cosmic variance on scales 
of one square degree (the size of our cones). We further point 
out that the mean LBG count derived from running a large 
number of realisation contains the full information available 
in the simulation. 

We ran 20 different light-cones, and found a mean num- 
ber of LBGs < NlbGs >= 4133 per square degree, with a 
standard deviation of o~n LBGs = 183. The modelled density 
of LBGs is thus 1.15 arcmin -2 , cons istent with the observed 
densit y of 1.22±0.18 arcmin -2 (e.g. IGiavalisco fc Dickinson! 
(2001)). In Fig. 21 we show the mean number counts, ob- 
tained from our 20 mock catalogues (solid line), as well as 
the dispersion, given at 1-cr by th e grey area. Our counts 
compa re nicely with the data from IGiavalisco fc Dickinson! 
(2001), over-plotted as diamonds. Note also that, as in real 
observations, we find very few interlopers (about 0.1%) with 
z < 2.3. The dashed line shows counts for our modelled 
LBGs when we do not apply any maturity selection, and fo- 
cus on a single example cone. The difference between these 
counts and the solid line shows that the history resolution 
mentioned in Sec . 12 . 2 . 21 mainly and strongly affects the faint 
end of the luminosity function (LF). Our most robust results 
thus concern the bright end of the modelled LBGs, where 
resolution does not affect galaxy properties. The dot-dashed 
line shows the counts for the LBGs which are older than 
1.3 Gyr. The dotted line shows the counts obtained from 
a mock observing cone of 0.09 square degrees built with a 
higher-resolution ACDM simulation. This higher-resolution 
simulation contains 320 3 particles in a cubic volume of side 
32/i _1 Mpc. It has identical cosmological parameters as given 



in Sec. 12.1.11 and the GALICS post-processing was done us- 
ing the same astrophysical parameter set as in GALICS I. 
The smaller volume and larger number of particles result in 
a galaxy mass resolution M res — 3.77 x 10 s M@, i.e. ~ 60 
times better resolution than in the fiducial simulation dis- 
cussed here. As expected, the introduction of a maturity age 
tmat is irrelevant for the LBG counts of this high-resolution 
simulation. The faint-end normalisation of the counts sup- 
ports our choice of tmat = 1-1 Gyr for our maturity crite- 
rion. The discrepancy in the bright counts are understood in 
terms of volume limit of the high-resolution simulation : just 
from Poisson statistics we expect about 30 times less objects 
per bin of luminosity since the volume is reduced by such a 
factor compared to the fiducial simulation. Furthermore, we 
have to account for cosmic variance which will affect more 
dramatically objects that are more clustered (see Sec. 13. 51 . 
Therefore, we cannot consider our high-resolution simula- 
tion as representative of the Universe as a whole, and for 
this reason, we prefer to use it only to fix the maturity cri- 
terion in the larger, low-resolution simulation, and we take 
this latter to make extensive predictions. 

In the remainder of the paper, except where mentioned, 
we will use a single realisation of a mock catalogue. Since 
we have no information on how the data might be biased 
(cosmic variance), we can choose any of our available cones. 
We settle for one which contains NlbGs = 4156, in good 
agreement with the data, and close to our mean value. 

3.3 Efficiency of the selection 

Internal properties of galaxies are not fully eclipsed by the 
IGM attenuation, so one cannot expect to extract all the 
galaxies at a given redshift, and only these, from a simple 
region selection in the colour-colour plane. It is therefore 
interesting to define an efficiency of the selection process, 
which should account for the fraction of lost galaxies - that 
is, observable galaxies with the right redshift which do not 
meet the colour criteria - and the fraction of interlopers 
restricted here to galaxies lying outside the redshift range 
[2.7; 3.4]. We therefore define the two following quantities: 

• the completeness (C), which is the ratio of the num- 
ber of selected galaxies with redshift ~ 3 to the number of 
detectable galaxies (i.e. galaxies with R < 25.5) at z ~ 3. 

• the confirmation rate (CR) which is the ratio of the 
number of selected galaxies with z ~ 3 to the total number 
of selected galaxies (with any redshift). 

With these two complementary quantities, we can define 
a measurement of the efficiency of a selection process as 
E = C x CR. This quantity takes values from to 1, mean- 
ing no source was detected at the expected redshift, and 1 
meaning that all the detectable sources with z £ [2.7; 3.4] 
were selected, and only these. For our sample of galaxies ex- 
tracted from a 1 square degree field, we find the following 
values for the colour criteria given by Steidel and co-workers 
(eqs|TJ : C = 96%, CR = 88%, and E ~ 85%. This high ef- 
ficiency shows that the LBG selection is indeed a very good 
method for selecting distant galaxies, and the high complete- 
ness shows that this selection grabs 96% of the detectable 
galaxies at z ~ 3. 

In an attempt to mimic observations, we also added 
photometric errors to our magnitudes. These are simply 
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modelled here as uniformly-distributed random values, be- 
tween ±<5 mag, in the three bands U, G, R. The selection 
efficiency drops to E — 70% when 8 reaches 0.25 mags. For 
a more conservative error amplitude of 0.1 magnitude, the 
efficiency settles around 80%. Obviously, photometric errors 
are far more complex than this simple model and should 
take into account systematics. Nevertheless, their impact on 
LBG selection as assessed by our simple estimate constitutes 
a useful first estimate. 

3.4 New contour 

Until now, the colour selection used to select distant galax- 
ies from multi-band imaging observations has been defined 
by placing synthetic galaxies, with a 'wide range' of ages, 
star forma tion histories, meta ll icities, dust content and red- 
shifts (e.g. lSteidel et all (ll995MMadau et al.Nl996l) '). onto a 
colour-colour plane, and by finding out which region encloses 
most distant galaxies and fewest interlopers. This method, 
although spectroscopic follow-up showed it gives good re- 
sults, sweeps under the carpet the question of the nature 
of the distant objects. It is expected to be robust anyway 
since IGM extinction -which is external to galaxies- plays 
a dominant role in the location of galaxies in the colour- 
colour plane. However, one might wonder whether a tighter 
set of criteria could be found if the 'wide range' of properties 
was replaced by more 'plausible' properties for galaxies at 
a given redshift. We present here the first attempt to build 
such criteria with an ab initio model of galaxy formation, in 
which distant galaxies' properties come out naturally of the 
hierarchical evolution within the dark matter content of the 
universe. 

We found our best contour in the colour-colour plane 
using the downhill simplex method jPress et alJll992T) in or- 
der to maximise the efficiency of the LBG-selection for a 
pentagonal contour with two sides bound to be vertical, and 
one horizontal in the (U n — G; G — R) plane. For this max- 
imisation, we used our full fiducial catalogue of one square 
degree so as to have a significant number of LBGs. The con- 
tour we calculate is the following : 
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This contour (the dashed lines of Fig. |3J gives an effi- 
ciency E — 88%, instead of E = 85% for the contour defined 
by Eq0 Our efficiency is higher because we discard a large 
fraction of high-z interlopers by setting an upper limit on 
the U n — G colour, which increases our confirmation rate to 
CR = 95%. However, our contour decreases the complete- 
ness (C = 94%), precisely because of this more accurate 
selection. In practise, the aim of LBG surveys is twofold : 

(i) select a large number of high-z candidates, with redshifts 
as high as possible, allowing for a broad redshift range, and 

(ii) use the detected galaxies to describe the properties of 
'normal' galaxies at a given redshift. In case (i), one wishes 
to favour completeness relative to confirmation rate. In the 
second case, one wishes to reach a compromise between C 
and CR. Bearing this in mind, it is worth pointing out that 



our selection is quite similar to that used by Steidel and col- 
laborators. This gives the overall impression that the UV- 
dropout technique is quite robust, and independent of stellar 
population modelling. As mentioned earlier, this robustness 
was expected since the locus of galaxies in the colour-colour 
plane is mainly given by IGM absorption, and thus by red- 
shift. 

3.5 Cosmic variance 

The degree to which observations of small areas are repre- 
sentative of the whole Universe is a crucial question when 
one wishes to interpret data in the cosmological context. 
Since our model contains spatial information, it is possible 
to study the cosmic variance of LBG densities, at least to 
a certain extent, because our simulation is volume-limited 
and does not include long-wavelength fluctuations. 

We estimated the cosmic variance as the field-to-field 
variation in the number of LBGs in sub-fields, extracted 
from our 20 different 1 deg 2 fields. For each size of sub-fields, 
ranging from 4 arcmin 2 to 1600 arcmin 2 , we used 20000 
random sub-catalogues (1000 per 1 deg 2 mock catalogue) 
to estimate the mean number of LBGs and the associated 
variance. 

In table^ we give the mean numbers of LBGs and the 
corresponding standard deviations for fields of 4, 16, 64,400 
and 1600 arcmin 2 . Numbers are given as well for LBGs 
brighter than R = 25 and R = 24.5. The fact that the 
standard deviation is higher by a factor > 2 to the pois- 
sonian deviation is a direct signature of the high clustering 
properties of LBGs. 

In Fig. |S] we show the variation of the relative devi- 
ation (a/{N)) versus the surface of the mock survey. The 
continuous line shows the measured deviation and the dot- 
ted line shows the Poissonian prediction. The dashed line 
shows the quadratic difference between those : <r c / (N) = 
(o-%~o-p oisson ) 1 '' 2 /{N}, where <r c is the contribution of clus- 
tering to the cosmic variance. The three crosses at 3600 
arcmin 2 show the total dispersion and the contributions of 
clustering and Poisson, from top to bottom. These were es- 
timated only from our 20 mock catalogues and thus rely on 
less statistics than the curves. On small scales (< 7 arcmin 2 ), 
where the sources are rare enough, the Poissonian noise dom- 
inates the variance. On larger scales, the variance is domi- 
nated by the clustering of LBGs, as the Poissonian variance 
falls off more rapidly. As the Universe is homogeneous on 
large scales, we expect the clustering contribution to dwin- 
dle progressively as the size of the survey increases. This 
regime is however out of reach with our current simulation 
because the simulated volume has an angular size of ~ 1 
square degree at z ~ 3. Above this scale we are bound to 
under-estimate the dispersion. Because of this, we cannot 
tell at the moment whether the points at 3600 arcmin 2 in 
Fig. |2| are low due to finite volume effect or to the fact that 
we reached the point where clustering strength starts to van- 
ish. 



4 CLUSTERING 

The clustering properties of LBGs give important con- 
straints on the dynamical processes that drive galaxy for- 



8 Blaizot et al. 



LBGs 


4 arcmin 2 
< N > (a N ) 
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< N > (<tjv) 


64 arcmin 2 
< N > (a N ) 


400 arcmin 2 
< N > (a N ) 


1600 arcmin 2 

< N > (<7 N ) 


R < 25.5 
R < 25. 
R < 24.5 


4.62 (2.76) 
2.71 (1.94) 
1.21 (1.20) 


18.35 (7.07) 
10.74 (4.80) 
4.78 (2.74) 


73.00 (18.9) 
42.74 (12.2) 
19.04 (6.69) 


454.6 (66.9) 
266.1 (42.1) 
118.3 (21.4) 


1817 (154) 
1064 (93.9) 
473.1 (46.6) 



Table 1. Mean expected numbers of LBGs brighter than R = 25.5, 25 and 24.5, in fields of 4, 16, 64, 400 and 1600 arcmin 2 . The standard 
deviations are also given between brackets. These standard deviations are systematically higher than the square-root of the number of 
sources, because of the strong effect of clustering. 
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Figure 5. Standard deviation in number counts of LBGs as a 
function of the survey area in square arcminutes. The solid curve 
is the measured standard deviation, and the dotted one is the 
Poissonian prediction. The dashed curve shows the quadratic dif- 
ference between the measured deviation and the Poissonian one, 
namely, the contribution of clustering to the cosmic variance. 

mation. Moreover, they can be observed directly and no as- 
sumptions are necessary to compute the angular correlation 
function. In this section, we show different aspects of the 
clustering properties of our modelled LBGs, namely, halo 
masses, the location of LBGs in the cosmological simula- 
tion, the angular correlation function, the spatial correlation 
function, and the bias of LBGs relative to dark matter. 



Figure 6. Masses of LBG host haloes as a function of redshift. 
The horizontal dotted line shows the median mass value of ~ 16 
10 U M Q . 
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Table 2. Median masses of DM haloes harbouring LBGs. 



4.1 Halo masses 

A first impression of the clustering properties of LBGs can 
be given by the mass distribution of the dark matter haloes 
which harbour them. In Fig. |S] we show the masses of these 
haloes versus the redshifts at which the LBGs are detected 
in our simulation. The median halo mass for our LBGs is 
15.8 x 1O 11 M0, as indicated by the horizontal dotted line. 
This mass is about ten times the halo mass resolution of 
the simulation. The median values for populations of mod- 
elled LBGs with various luminosities are given in table |5] 
The relatively large mass of the haloes harbouring our mod- 
elled LBGs at 2 ~ 3 indicates that these galaxies are highly 
biased with respect to the dark matter distribution. This 
matches qualitatively what is observed in diff erent surveys 
jGiavalisco et alJ <1998h . lArnouts et alJ tOO'j) ). 



4.2 Location of LBGs and their descendants 

On the left-hand side panel of Fig. [7] we show the location of 
LBGs at z ~ 3 in a slice of our simulated cosmological vol- 
ume. The white discs show LBGs, and the underlying grey 
background shows the dark matter density (the darker the 
denser), in logarithmic scale. The LBGs are clearly highly 
clustered in this picture. The right-hand panel shows the 
location of the descendants of LBGs at z — as squares. 
Black squares indicating spiral galaxies, and white squares 
elliptical or lenticular galaxies. The white discs are ellipti- 
cal/lenticular galaxies which do not have a LBG progenitor 
at z ~ 3. The descendants of LBGs are even more highly 
clustered than LBGs themselves, as they fell into DM struc- 
tures from z = 3 to z = 0. We will come back to the prop- 
erties of the descendants of LBGs in section ItT^l 
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Figure 7. Le/i /iartd panel : location of LBGs at z=3. White discs show LBGs, and the DM density (in log) is shown in grey scale, the 
darker the denser. Right hand panel : location of local ellipticals and the descendants of LBGs in the same comoving volume, at redshift 
0. The white symbols show elliptical or lenticular galaxies, the squares indicate descendants of LBGs that are present in the volume of 
the left hand side panel. Black squares show spiral descendants of LBGs. The LBGs and their descendants are highly clustered, located 
in dense DM regions. The projected volume shown is 150 Mpc on a side, and has a depth of 15 Mpc. 



4.3 Angular correlation function 

For the computation of the angular cor relation function 
(ACF ), we used the estimator proposed bv lLandv fc Szalavl 
fl993|) (hereafter LS93) : 



w(8) 



DD{0) - 2DR(8) + RR(8) 
RR(6) ' 



(2) 



where DD(6) is the number of pairs of LBGs with angular 
separation between 8 and 8 + 88, RR(0) is the analogous 
quantity for a random catalogue, and DR(8) the number of 
observed-random pairs. 

The estimate of the errors in the general case is a dif- 
ficult task. Here, we consider Poissonian noise (in the num- 
ber of pairs) only. We expect this to be a good approxima- 
tion because our survey is large enough (compared to the 
scales we probe), and we deal with a rather diluted distribu- 
tion of galaxies so that finite- volume errors will be negligible 
compared to discreteness errors. Moreover, considering only 
Poissonian error-bars implies that we underestimate the real 
uncertainties. It is not a concern here as our goal is to show 
that our spatial distribution of LBGs is consistent with the 
observed one - and this is the case within Poissonian uncer- 
tainties (see next section). Following LS93, we write Poisso- 
nian errors as : 

\1 1/2 



l+w{8) 



1) 



N r (N r 



RR 



(3) 



where N g is the number of galaxies in our mock survey, and 
N r the number of objects in our random catalogue (here, 
N r ~7x 10 4 ). The first factor in the above equation shows 
the Poissonian contribution, as it is inversely proportional 
to the square root of the number of LBG pairs. The sec- 



ond term is a geometrical contribution, and depends only 
on the shape of the survey. In this expression of the stan- 
dard deviation, we have assumed that the integral constraint 
was negligible. Of course, the standard deviation of Eq. 
strongly depends on the angular bin-size we use to evaluate 
the ACF. This appears only implicitly : if the bin size dou- 
bles, RR will also roughly double, and so a w will be divided 
by a factor %/2. I n ord er to compare our results to those of 
iGiavalisco et alj il99Sft . we adopt the same binning as these 
authors. 

Finally, note that the random transverse shifting of 
boxes every 150 comoving Mpc along the line of sight that we 
use to suppress replication effects implies an underestimate 
of the ACF of about 10% at most, as detailed in momaf. 

Our measured ACF is represented in Fig. |H| by crosses 
and their error-bars. The solid line shows a maximum likely- 
hood fit of a power-law w(9) = A w 8~^ to our data. The best 
fit parameters are : A w = 4.41±0.7 8 and = 0.9±0.04. Thi s 
is consistent with the results from iGiavalisco et alJ (l998). 
Their fit to observed data is shown in Fig. |S]by the shaded 
region 1 . The vertical dotted line roughly shows the angular 
size of a cluster of galaxies at z ~ 3, below which the con- 
tribution to clustering mainly comes from pairs of galaxies 
within the same halos. In GALICS, the positions of galaxies 
inside halos are given according to some prescription and 
not directly measured from the N-body simulation. Hence 
our model only makes reliable predictions for the two-halos 
contribution to w (8). 



1 The values from these authors that we consider are those given 
by their so-called "LS-random" estimate, which corresponds to 
our estimate. 
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Figure 8. Angular correlation function for our sample of LBGs, 
extracted from a 1 square degree mock catalogue. The crosses 
show our estimate of w(8) with the LS93 estimator, and the error 
bars come from a Poissonian estimate. The grey area shows data 
from lGiavalisco et alj ll99Sl) ■ This area corresponds to their "LS- 
random" estimate, which is similar to ours. The solid line shows a 
maximum-likelihood fit to our data by a function w(9) = A w 9~ P, 
with best fit values indicated in the panel. 



4.4 Spatial correlation function and linear bias 

We computed the spatial correlation function (SCF) using 
the same estimator (LS93) as for the ACF. This time, how- 
ever, we computed the SCF directly from the snapshots of 
GALICS, bypassing the cone of sight generation. We show the 
result in Fig. El The top panel is the SCF of LBGs identified 
in the snapshot at z = 3. The error bars only show the pois- 
sonian deviations. The bottom panel shows the linear bias of 
LBGs relative to the dark matter particles of the cosmologi- 
cal N-body simulation (6 2 (r) = £e,bg(t") /£dm(t)). We limit 
our plots to the approximate range 1-10 h~ Mpc for two 
reasons : (i) we distribute galaxies inside clusters accord- 
ing to a semi-analytic spherically symmetric prescription, 
so we loose angular information for some galaxies below 1 
h~ Mpc, and (ii) to avoid edge effects, we limit the study 
to scales below one tenth of the simulated volume size (100 
ft -1 Mpc). 

We find that the linear bias slowly decreases from ~ 5 
to ~ 3 between 0.5 and 20 comoving h^Mpc. This is com- 
patible with the bias es tima ted bv iGiavalisco et al.l 11998), 
lAdelberger et all (Il998l) and lPorciani fc Giavaliscall2002lL 



5 UV/IR LUMINOSITY BUDGET 

One of the most important challenges toward understand- 
ing galaxy formation is to determine at which epoch stars 
formed, or, in other words, how and when did the mass of 
stars we see in local Universe galaxies assemble. The ob- 
servation of distant galaxies in the optical window is now 
known to be insufficien t to retrieve star formation rates (e.g. 
iDevriendt et alJ lll999ll 'l. The discovery of the cosmic far in- 
frared background, and the numerous IR/sub-mm surveys, 
have made it clear that dust extinction within galaxies plays 
a fundamental role in determining the UV emission of dis- 
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Figure 9. Top panel : spatial correlation function of our modelled 
LBGs in the snapshot corresponding to z=3. Bottom panel : linear 
bias of our modelled LBGs with respect to the dark matter. 



tant galaxies, and therefore in extracting information on the 
cosmic star formation rate. The reason is that the UV lu- 
minosity of a galaxy is governed by the amount of recent 
star formation -which is responsible for young massive stars 
dominating the energetic part of the SED- and dust ab- 
sorption -which is most efficient at UV wavelengths, and 
re-emits light in the IR/sub-mm domain. Unfortunately, far 
infrared or sub-mm instruments do not yet have sensitiv- 
ity or resolution which would make a thorough follow-up of 
LBGs at large wavelengths possible. 

Our model includes a treatment for dust extinction and 
emission in the IR/sub-mm domains, which allows us to pre- 
dict extinction properties of LBGs, and make the connection 
between LBGs, seen in the optical, and sources detected 
at IR/sub-mm wavelengths. In this section, we first discuss 
the UV luminosity function and selection effects. Second, 
we consider the extinction properties of these objects as 
predicted by GALICS. Third, we s how how both effec ts al- 
ter the so-called Madau diagram jMadau et al.l ll996T We 
eventually make predictions concerning the detectability of 
the sub-mm counterparts of LBGs. 



5.1 Luminosity functions at z ~ 3 

In Fig. 1101 we compare different estimates of our measured 
luminosity function (LF) at z ~ 3 with the LF given by 
ISteidel et all ll 19991) . 



• The solid histogram shows t he LF for all our mo delled 
LBGs, computed in the same wav lSteidel et all 1^999) com- 
puted their observed LF. First we divide the number of 
LBGs in each apparent-magnitude bin by the effective vol- 
umes quoted by these authors for our cosmological param- 
eters. Then, we convert our apparent magnitudes to abso- 
lute magnitudes using the relation L v = iq- -< as - & +™ab) x 
47rd£/(l + z) where d L = 5.6 10 28 h' 1 cm, is the s i ngle v alue 



of the luminosity distance given bv lSteid'el et"ail i 19991) for 
2 = 3.04 in a ACDM universe (Q A = 0.7, and VL m = 0.3). 
Note that the difference in luminosities induced by using a 
single redshift and luminosity distance are negligible when 
compared to the uncertainties in the volume corrections. 
• The dashed histogram shows the luminosity function of 
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Figure 10. Luminosity functions at z ~ 3, for LBGs. The x-axis 
shows absolute R magnitudes in the observer fra me. The black 
diamonds show the data from ISteidel et all Jl999l) . converted to 
the current cosmology. The solid histogram shows the LF we com- 
pute for our sample of LBGs, using the effective volumes given in 
ISteidel et alJ il99fll) . We assume a single luminosity distance for 
the whole sample (also given by ISteidel eta?] <1999ft 1. This his- 
togram is in fair agreement with the data. The dashed histogram 
shows the LF computed for all the detectable galaxies (i.e. galax- 
ies with R < 25.5) in our sample that have a redshift between 
2.7 and 3.4. Here, we normalised the LF to the total volume of 
the cone of sight between redshifts 2.7 and 3.4. We also used the 
"true" luminosity distance of each galaxy. This LF is lower than 
the previous one because of volume corrections. Finally the dotted 
histogram shows the true luminosity function of LBGs only, i.e. 
without volume corrections. The two vertical dotted lines show 
the resolution and detection limits. 



all detectable galaxies (R < 25.5) with redshifts between 2.7 
and 3.4. This time, we first convert apparent magnitudes to 
absolute magnitudes using each galaxy's luminosity distance 
and redshift. We then divide the number of objects in each 
magnitude bin by the total volume of the cone of sight be- 
tween 2 = 2.7 and z = 3.4, ignoring any volume effect due 
to the selection function. 

• The dotted histogram is the same as the previous one, 
but only includes LBGs. It lies just below the dashed his- 
togram in every bin, showing that LBG selection is more or 
less equally efficient at all magnitudes. 

The resolution and detection limits are shown as dotted ver- 
tical lines. The f ormer limit was defi ned in section The 
latter is given bv lSteidel et alJ lll999l) and corresponds to the 
apparent magnitude limit R < 25.5. Our resolution limit is 
of the same order. 

The solid curve of Fig. llUI which reproduces the analysis 
of ISteidel et alJ <ll999l) . yields a very encouraging match to 
the data. At this stage, we remind the reader that in this 
paper we use the same galics model that was used in 
GALICS I to reproduce local galaxy properties. It is therefore 
an impressive result that we manage to match the observed 
luminosity function (or the number counts shown in Fig. 
of galaxies at high redshift as well. 

In Fig. 1111 we show the effect of dust extinction on lu- 
minosity functions. The solid histogram on this figure repre- 



Figure 11. Effect of dust extinction on the luminosity function 
of galaxies at z ~ 3. The solid histogram shows the LF of all 
our modelled LBGs at z ~ 3, as a function of their apparent R 
magnitudes. The dashed histogram is the LF these galaxies would 
have if they contained no dust. It was computed directly from the 
rest-frame emission of galaxies at 1600A, assuming no extinction. 
Finally the dotted histogram presents the same results as the 
dashed one but shifted by 2 magnitudes to facilitate comparison 
with the solid histogram. 



sents the absolute R magnitude (in the observer frame) lu- 
minosity function, and includes all our modelled galaxies at 
z ~ 3 (with the usual R < 25.5 criterion on apparent magni- 
tude). The dashed histogram corresponds to the luminosity 
function of all these galaxies computed directly from their 
luminosities at 1600A in the rest-frame, assuming no extinc- 
tion. The dotted histogram is the dashed one shifted two 
magnitudes fainter to facilitate the comparison with the ex- 
tinguished LF (solid histogram). Note that the filters differ 
for the two LFs. Indeed, the R filter blue-shifted to z ~ 3 has 
a width AA ~ 375A and is centred at A_r ~ 1700A, whereas 
the filter we used to compute the 1600A emission is centred 
at this wavelength (in the galaxy rest-frame), and is a top- 
hat function of width 20 A. Nevertheless the huge difference 
between the two histograms demonstrates how delicate it is 
to extract the UV flux really emitted from observed broad 
band fluxes. We emphasise that it is this very same UV flux 
which is mainly used to compute the cosmic star formation 
rate at intermediate and high redshift. Moreover, as a final 
note of caution, we stress that even the slopes of the LFs 
differ, which means that a unique conversion factor in not 
sufficient to accurately correct data for extinction. 



5.2 Extinction of LBGs 

Obscuration by dust plays a key role for the extraction of 
star formation rates from UV observations. Here, we show 
how LBGs and other detectable galaxies at 2 ~ 3 are af- 
fected by this process, and we give the mean extinction fac- 
tors predicted by GALICS for these galaxies. 
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Figure 12. Ratio of IR luminosity (8^m < A < 1000 fim) to 
bolometric luminosity for LBGs (dots) and detectable galaxies 
with 2 G [2.7; 3.4] (crosses) that do not meet the LBG selection 
criteria. Most of these galaxies emit more than half their light in 
the IR/sub-mm. 



5.2.1 



Qualitative view 



In Fig. 1121 we show an estimate of the internal extinction 
versus redshift for our sample of LBGs and other galaxies 
at z ~ 3 which are brighter than R — 25.5 (the so-called 
'detectable' galaxies). Internal extinction is estimated here 
as the ratio of infrared bolometric luminosity (integrated 
from 8/im to 1000/im at rest) to bolometric luminosity. A 
value of 0.5 of this ratio means that half the light emitted 
by stars is absorbed and re-emitted by dust. We find that 
~95% of our LBGs emit more than half their light in the 
IR/sub-mm. This fraction of heavily extinguished objects 
holds for all detectable galaxies at z ~ 3. This suggests that 
an important contribution of the IR/submm background is 
due to z ~ 3 galaxies. 

The fact that missed galaxies (crosses in Fig. 1121 lie 
on rather narrow redshift bands on either side of z ~ 3 is 
another manifestation of the LBG selection efficiency: fore- 
ground galaxies are only missed because of border effects 
in the colour-colour plane, whereas background galaxies are 
more extinguished on average. 

5.2.2 Extinction factor 

The usual way to quantify extinction is through the absorp- 
tion coefficient A{\) = k(X)E(B — V), where k(X) is the 
adopted attenuation law. In terms of luminosities, this coef- 
ficient can be written as 

Lno ext\ 



A(X) = 2.5 log - 



t(A) 1 



(4) 



where L no ex t(X) is the non-extinguished luminosity at 
wavelength A, and L ext (X) the extinguished one, computed 
for a given inclination of the galaxy. The extinction fac- 
tor is then defined by F(X) = 10 4A(A) , that is, in terms 
of luminosities, F(X) = L no ext(X)/ L ext (X). This quantity 
is of great importance because it is necessary for the ex- 
traction of information about the star formation rates of 



LBGs 


F e// (1600) 


F a „ e (1600) 


R < 25.5 


6.19 


4.75 


R < 25 


6.13 


4.67 


R < 24.5 


6.21 


4.83 









Table 3. Extinction factors at 1600 A for different sub- 
populations of modelled LBGs. The two estimators described by 
Massarotti et al (2001) are given for each sub-sample. 



distant galaxies. However, its value is still a matter of de- 
ba te. This is partly due to averaging effects, as emphasised 
bv lMassarotti et alJ <l200ll) . Indeed these authors stress that 
computing an average extinction factor from an average ab- 
sorption coefficient - namely 



F ave {X) = 10 



oaJ2- A iW/ N 



10 



0.4<A(A)> 



(5) 



where N is the number of galaxies in the sample - leads to 
an underestimate of the effect of dust. They suggest that 
one should rather use an effective extinction factor defined 
by 



Fe//(A ) = W 0AA ' W = 10°- 4A ^ (A) . 



(6) 



In order to understand the difference between both esti- 
mates, it is easier to re-write them in terms of luminosities. 
The average extinction factor becomes the geometric mean 

- l/N 



Fave(X) — 



n 



*-><r 



whereas the effective extinction factor is simply the arith- 
metic mean 



t (A) 



t (A) 



(8) 



We have computed both factors for our modelled galaxies, in 
order to facilitate the comparison with lSteidel et al] (Il999l) 
who use _F a „ e (1600) = 4.7, and also to see whether they 
give us results which d iffer by a factor ~ 3, as claimed by 
iMassarotti et all (1200 ll) . We summarise our results in table 
13 We only find a discrepancy of about a factor 1.3 between 
the two estimators, with extinction factors F e /f ~ 6.2 and 
Fave ~ 4.75. These values are much lower than the value 
of 12±2 obtained by Massarroti and co-workers, but their 
analysis is based on the Hubble Deep Field only, and in- 
cludes galaxies spanning a much wider range of redshifts 
than those we consider here. Our average extinction factor 
has a value close to that used bv lSteidel et al] (^999), and 
our effectiv e extinction factor is c o nsiste nt with the detailed 
analysis of lAdelberger fc Steidell ]2000l) who conclude the 
effective extinction should be around a factor ~ 8. 



5.3 Cosmic star formation rate 

Another uncertain step in retrieving star formation rates 
from UV fluxes is the UV-SFR connection, considering no 
absorption (or no dust). This relation is quite sensitive to 
the initial mass function (IMF) and we remind the reader 
that the IMF in this analysis is the same we used in our 
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LBGs 


SFRioo/ L(1600) 


SFR inst / L(1600) 


R < 25.5 


16.7 (3.3) 


8.94 (2.2) 


R < 25 


16.8 (3.6) 


8.88 (1.1) 


R < 24.5 


16.8 (4.5) 


8.77 (0.9) 



Table 4. SFR to UV luminosity ratios (in Mq yr _1 per 10 29 erg 
s _1 Hz -1 ) computed for different subsamples of our modelled 
LBGs. The UV luminosities are computed as if the galaxies con- 
tained no dust, that is, with no extinction. In the first column, 
the SFR is computed as the mass of stars younger than 100 Myr 
divided by 100 Myr. In the second column, the SFR is the in- 
stantaneous star formation rate, computed over a Myr only. For 
sake of completeness we also give standard deviations between 
brackets for each sub-class of LBGs. 

fiducial model described in GALICS I, namely a iKennicuttl 
fllQSot IMF from 0.1 M Q to 120 M©. 

5.3.1 UV/SFR conversion 

The first step is to obtain extinction factors, in order to re- 
late the observed magnitudes to the UV flux that galaxies 
emit prior to extinction. The second step is to make the link 
between these extinction-corrected UV luminosities and star 
formation rates. This effectively is feasible because massive 
stars are the main contributors to the galaxy UV flux and 
are sufficiently short-lived. However, star formation rates 
have to be averaged over a period of about 100 Myr so that 
a robu st enough relation may be establish (e.g. IKennicuttl 
(1998)). In tableHwe give SFR-to-UV luminosity ratios for 
different sub-populations of LBGs, whose SFRs are evalu- 
ated as the mass of stars younger than 100 Myr, divided by 
that period of time. Because of the rapid timescale of stel- 
lar evolution, this is an under-estimate of the real SFR, by 
about 20% for Kennicutt IMF. We also give ratios computed 
using the instantaneous SFR for comparison. Both these es- 
timates yield reasonable values because the evolution of the 
SFR is generally sm all over a 100 Myr p eriod and therefore 
are consistent with iMadau et al. ( 1996) - who use a con- 
version factor of 9.54 Mq yr _1 per 10 29 erg s _1 Hz -1 - and 
with Kennicutt (1998), who gives a value of 14 Mq yr _1 
per 10 29 erg s _1 Hz -1 . The ratios we compute, however, are 
the result of particular star formation histories: those of our 
LBGs. Looking at table^] one remarks that averaged SFRs 
are systematically larger than their instantaneous counter- 
parts by about a factor two. This indicates that modelled 
LBGs tend to be identified as such when they enter a "post- 
starburst" phase. 

5.3.2 The cosmic star formation rate history 

Following ISteidel et alJ l^!99 : ), we compute the cosmic SFR 
at z ~ 3 from our UV luminosity function. Integrating the 
modelled UV LF of LBGs and converting it into a star for- 
mation rate density using Table0] we find psfr = 8.9 x 10 -3 
M0yr -1 Mpc -3 . Correcting this value for extinction with the 
effective factor of Table[3]yields a density psfr = 5.5 x 10 -2 
M0yr _1 Mpc -3 . The last step is to correct for the missing 
faint end of the luminosity fu nction (unresolved galaxies). 
Following ISteidel et alJ l)l999t) . we apply a correction of 0.5 



dex and fin d psfr = 0- 1 7 M^ yr * Mpc 3 , which is con- 
sistent with ISteidel et al.l lll999l) and lAdelberger fc Steidell 
J2000D . 

5.4 Prediction of sub-mm flux 

One of the main ambitions and major difficulties of hierar- 
chical models of galaxy formation is to reproduce the sub- 
mm number counts. Several attempts have been made in 
semi-analytic models to fit these observations, but none of 
the recipes has yet em erged as a well accepted solution. 
iGuiderdoni et all dl998h implemented a massive IMF to ob- 
tain SCUBA objects, but there is neither a real physical 
motivation for this solution, no r strong observational hints. 
iDevriendt fc Guiderdonil j2000l) introduced an ad-hoc red- 
shift evolution which did n ot possess firmer foundations. Fi- 
nally, iBalland et al.l i2003T) used a simple yet physically mo- 
tivated model to estimate energy exchanges during galaxy 
interactions which they propose as the principal mechanism 
for triggering starbursts. They show that this solution yields 
a good match to mid-IR as well as SCUBA counts. A thor- 
ough discussion on this topic will be given in paper 4, but 
we forewarn the reader that the fiducial GALICS model pre- 
sented in GALICS I that we are using for the present study 
probably underestimates the cosmic amount of 850 pm light 
by a factor ~2. 

In Fig. ESI we show ratios of IR to UV (1600A) lumi- 
nosities versus the sum of these luminosities for our mod- 
elled LBGs. As a crude approximation, this is equivalent to 
"dust absorption" versus "star formation rate" . The grey re- 
gions broadly indicat e the c ontour of the data analysed by 
lAdelberger fc Steidell i2000t) . and the solid line shows the 
SCUBA detection limit. The good match between obser- 
vations and our modelled LBGs suggests that we manage 
to compute dust absorption and star formation rates in a 
consistent way. We also recover the observed correlation be- 
tween star formation and dust content of the galaxies : the 
more obscured, the more star-forming. However, we predict 
that only about 1 % of our modelled LBGs should be de- 
tectable with SCUBA at 85 pm at the 2 mjy level , in agree- 
ment with the estimates of lChapman et all (|2000T) . 

The model can be used to make predictions for other 
IR/submm surveys. For instance, we find that the SWIRE 2 
survey with SIRTF will be able to detect only 0.04% (resp. 
0%, 0%, 0%, 0.02%, 1.4%, 2.3%) of them at 170 pm (resp. 
70pm, 24pm, 8pm, 5.8pm,4.5pm, 3.6pm) at the 17.5 mjy 
level (resp. 2.75mJy, 0.45mJy, 32.5/^Jy, 27.5/^Jy, 9.7,/xJy 
7.3/iJy). These objects are clearly a key target for ALMA 
which should be able to detect all of them at 850 pm or 1.3 
mm, at the 0.1 mjy level. 

Predictions at other optical and IR/submm wavelengths 
can be easily retrieved from queries to our web interfaced 
relational database located at http://galics.iap.fr 



6 NATURE OF LBGS 

As emphasised bv lShaplev et all ll200ll) . the question of the 
nature of LBGs remains quite open. Different models sug- 

2 |http : //www . ipac ■ caltech.e du/SWIRE/ 
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Figure 13. Ratio of IR to 1600A luminosity versus the sum of 
these luminosities. In a first approximation, this shows absorp- 
tion by dust versus star formation rate. T he grey region ro ughly 
indic ates the location of observed LBGs jAdelberger &: Steidell 
2000), whereas the dots represent our modelled LBGs. The solid 
line marks the SCUBA detection limit at 850 jim. 



gest that LBGs are seen because they are undergoing strong 
starbursts due to major mergers, whereas others claim that 
LBGs are central galaxies of massive haloes that form stars 
steadily so as to reach a consequent size by redshift 3. The 
increasing number of observations seems to indicate that 
LBGs have many facets. Each of the two above-mentioned 
scenarios can explain some of their properties. We describe 
in this section the general properties of our modelled LBGs 
and investigate their star formation histories. 



6.1 Physical properties 

Sizes 

In Fig. 1141 we show the distribution of sizes of our modelled 
LBGs. The size of a galaxy is defined here as the stellar- 
mass weighted average of the three component (disc, bulge, 
starburst) half-mass radii. We find a median value of 2.6 kpc 
for our LBGs, which is in good agreement with the analysis 
bv lGiavalisco et "all il99fJl of a few LBGs observed with the 
HST (but identified with ground-based images) for which 
they find a typical size of about 2 kpc. The mean values and 
standard deviations are summarised in table |6~D 



Stellar masses 

The distribution of stellar masses of our modelled LBGs 
is given in Fig. 1141 The median stellar mass is 2.6 x 10 10 
Mq , in very good a greement with the me an stellar mass de- 
duced from data bv lShaplev et al] i200lT) (< M„tar >med= 
2.4 x 10 10 M©), and iPapovich et al] feOOlh (~ 10 10 M ). 
The shape of our distribution is consistent with that given 
by Shapley et al. (2001, their Fig. 10 b), and also with th e 
"accelerated quiescent" model of SPF dPrimack et all200ll) . 



Metallicity 

The distribution of metallicities of our modelled LBGs is 
given in Fie. 1141 The solid histogram shows the 3-component 
average ISM metallicity of our modelled galaxies, the dotted 
histogram that of the discs, the dashed line shows that of 
the bulges and the dot-dashed one that of the starbursts. We 
find a median metallicity of ab out Zp, / 5 , which is ag ain in 
good agreement with observations iPettini et alj EoOL). The 
starburst and bulge components have higher metallicities, 
due to the more intensive star formation taking place there. 

Star formation rates 

In Fig. 1141 we show the distribution of star formation rates 
(SFR) for our modelled LBGs computed in two different 
ways. The solid line shows the instantaneous star formation 
rate (SFRi„ st ), computed over the last Myr. This quantity 
is not directly comparable to the value deduced from obser- 
vations of_the_Uy throu gh broad band filters (see discus- 
sion in iKennicutd Jl998)). However, they indicate the cur- 
rent state of LBGs, and should compare better to estimates 
of SFR obtained from emission lines. For this estimate, we 
find a median value of 27.4 Mq yr _1 , lower than the ob- 
served value by a factor ~ 3. The dotted curve shows the 
distribution of star formation rates computed as the mass of 
stars younger than 100 Myr, divided by 100 Myr (SFR100). 
Th is is much closer to what is deduced from observations 
by IShaolev et al] l)200l]) as it takes into account the star 
formation history during a time of order the lifetime of B 
and A stars. Our measure is however an under-estimate of 
the real SFR (by about 20%) because we miss all the mass 
lost through stellar evolution during a 100 Myr. The median 
value of SFR100 is 50.8 M yr"\ with 20% of the modelled 
LBGs forming more that 100 Mq per year. These SFRs are 
in better agreement with the data, alth o ugh s till low com- 
pared to the estimate of IShaplev et al] 11200 1]) : ~ 90 Mq 

yr- 1 . 

The question of what triggers these h igh SFRs has 
been investigated by ISomerville et all i200ll) who suggest 
that star formation in LBGs is triggered by mergers. In our 
model, we find that only 30% of the LBGs have undergone 
a merger (although this figure goes up to 60% for our high 
resolution simulation). Thus we conclude that star forma- 
tion in LBGs is both due to merger-triggered starbursts and 
active cooling of gas on the galaxies. 

Conclusion 

The main point that the previous statistics make is that 
GALICS does not reproduce LBG counts by chance, as it 
reproduces their dynamical and chemical properties quite 
well. Again, we emphasise that this is a natural outcome of 
a model which also reproduces quite well the properties of 
local galaxies, as shown in GALICS I. 

6.2 The descendants of LBGs 

In a recent paper, iNaeamine] (1200211 investigated the nature 
of the progenitors and descendants of LBGs within a cos- 
mological hydrodynamical simulation. We extend here this 
investigation to the new methodology we have developed 
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Figure 14. Top left panel : Radius (kpc) distribution for our modelled LBGs. The radius is a stellar-mass weighted average radius of 
three components: disc, bulge and starburst. The median value of ~ 2.6 kpc, indicated with the dotted vertical line, is in good agreement 
with the data from the HDF iGiavalisco et ai]|l996h . Top right panel : Stellar mass distributions for our modelled LBGs. The median 
value of 2.6 X 1O 1O M0 is indicated by the dotted vertical line. Bottom left panel : The solid histogram shows the 3-component average 
ISM metallicity distribution of our modelled galaxies, the dotted histogram that of the discs, the dashed histogram that of the bulges 
and the dot-dashed one that of the starbursts. Bottom right panel : The solid histograms represent the distribution of instantaneous 
SFRs of our modelled LBGs, whereas the dotted histogram shows the distribution of SFRs averaged over 100 Myr for the same galaxies. 
The median SFRi nst is ~28 Mq yr _1 , as indicated by the left hand side vertical dotted line, and the median of SFRioo is ~50.8 Mq 
yr — 1 , indicated by the right hand side vertical dotted line. 



LBGs 


< ri /2 > a 
(kpc) 


< Mstar > b 

(10 10 M Q ) 


<Z> c 
(Z@) 


<SFRi nst >^ 
(M Q yr- 1 ) 


<SFRioo > e 
(M Q yr" 1 ) 


R < 25.5 


2.61 


2.60 


0.219 


27.4 


50.8 


R < 25 


2.61 


2.52 


0.221 


26.2 


48.4 


R < 24.5 


2.59 


2.20 


0.217 


23.8 


43.2 















a Half- mass radius of largest galaxy component (disc, bulge, or burst) (standard deviation). 

^ Stellar mass of galaxies (standard deviation). 

c Mean metallicity (standard deviation). 

^ Instantaneous star formation rate summed over the three components of our modelled galaxies. 

e Star formation rate summed over the three components of our modelled galaxies, and averaged over 100 Myr. 

Table 5. General properties of our modelled LBGs at z ~ 3. 



16 Blaizot et al. 



with GALICS and momaf. We first start with examples of 
star formation histories of modelled LBGs, following them 
until z = 0, and continue with an investigation of the nature 
of the local descendants of these galaxies. 



6.2.1 Star formation histories 

In Fig, 1151 we show several examples of star formation histo- 
ries for our modelled LBGs. The top panels show the star for- 
mation rate as a function of cosmic time; the middle panels 
the stellar (solid curves) and gas (dashed curve) mass evolu- 
tion; and the lower panels that of the metallicity. The verti- 
cal dotted lines indicate the moment at which each galaxy is 
identified as a LBG. Before this vertical line, the given his- 
tory is that of the most massive progenitor at each merger 
undergone by the future LBG. Each column represents a 
galaxy, and they are sorted from left to right in order of 
decreasing redshift of identification. 

We find a large variety of star formation histories, rep- 
resentative of the variety of states in which galaxies can be 
when they meet the LBG criteria. This also illustrates the 
two processes that trigger star formation in LBGs : active 
cooling of the hot gas in halos onto the young galaxies as 
opposed to merger-triggered starbursts. 

Finally, notice how some LBGs evolve quiescently until 
the present day, resulting in massive spirals of which they 
form the bulges, and how some undergo several extra merg- 
ers, resulting in present-day massive ellipticals. From the 
middle panels of Fig. 1151 it can be seen that our LBGs rep- 
resent, on average, about 10% of the total mass of their 
descendants in the local Universe. 



6.2.2 LBGs as progenitors of ellipticals ? 

We define the morphology of galaxies at z = using the 
ratio of bulge to disc B-band luminosities, as described in 
GALICS I. With this definition, we find that ~77% of our 
modelled LBGs end up in ellipticals or lenticular galaxies, 
the remainder ending up in spirals. On the other hand, we 
also computed the number of elliptical or lenticular galaxies 
at z — that have had at least one LBG progenitor in their 
merging history tree. These galaxies represent about 35% 
of the whole local population of early-type galaxies (with 
M B < -18). 

In Fig. 1161 we show the expected location of the de- 
scendants of our modelled LBGs. The x-axis is the size (i.e. 
number of galaxies) of groups or clusters of galaxies, and 
the y-axis gives the number of LBGs that end up in groups 
of these sizes. Elliptical or lenticular descendants of LBGs 
are found in groups with a median number of ~15 galaxies, 
while spiral descendants are found to be satellite galaxies in 
clusters/groups of ~30 galaxies. 



7 DISCUSSION 

7.1 Summary of results 

This paper is the third one of the GALICS series that presents 
a consistent attempt to model hierarchical galaxy formation 
within a hybrid approach in which dark matter collapse is 
computed in a large cosmological simulation, and the physics 
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Figure 16. Number of galaxies in haloes where LBG descendants 
are found, for three morphological types of descendants. Vertical 
dotted lines indicate median numbers. 



of baryons is described by semi-analytic recipes. The first 
paper (GALICS i) has detailed the model and the predicted 
statistical properties of local galaxies, which have been com- 
pared against data. The second paper (GALICS ii) has pre- 
sented the evolution of a set of these properties from red- 
shift 2 = 3 to z = within the hierarchical paradigm. In 
the current paper, we have proposed a more detailed study 
of the properties of optically-bright, star-forming galaxies 
at z ~ 3. We have converted the predictions of our GAL- 
ICS model into a mock observing cone by means of the MO- 
MAF package. From this observing cone, a sample of LBGs 
has been extracted with selection criteria as c lose as possi- 
ble to those used in the observational sample of lSteidel et alJ 
( 1996). This method allows us to get a three-fold insight on 
(i) the efficiency of the selection criteria to capture z ~ 3 
galaxies, (ii) the nature of the objects that have been se- 
lected according to these criteria, and (iii) the properties 
that could be measured in forthcoming follow-up. The re- 
sults of this paper have been obtained from the same cos- 
mological simulation (a ACDM model with 256 3 particles 
in a cubic box with lOO/i -1 Mpc on a side), and with the 
same astrophysical parameters in the semi-analytic post- 
processing as in the first two papers. In doing so, we are 
attempting to progressively build up a consistent view of 
hierarchical galaxy formation from z — 3 to z = 0. 

In this prospect, the successes of the model are satis- 
factory, in spite of several shortcomings. Our model LBGs 
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Figure 15. Examples of galaxy histories in a ACDM cosmology. Top panels : instantaneous SFR versus cosmic time. Middle panel 
: masses of stars (solid curves) and gas (dashed curves) in our LBGs versus cosmic time. Lower panels : evolution of average ISM 
metallicities with cosmic time. The vertical dotted lines indicate the time at which galaxies were identified as LBGs. 



were selected according to colour and magnitude criteria in 
the mock observing cone. This allowed us first to examine 
the efficiency of the selection method. We have seen that: 

• The model is able to reproduce the number density 
counts of LBGs on the sky, and gives a number density of 
-~1.2 arcmin -2 at the magnitude limit of the sample, in good 
agreement with the data. 

• The set of colour an d m agnitude criteria des igned by 
ISteidel fc Hamilton! lll993ft and lSteidel et alJ lll995l) to select 
z ~ 3 galaxies on the basis of synthetic stellar populations 
evolving according to the monolithic model, is well suited 
to the selection of LBGs within the hierarchical model. Our 
best contour in the colour-colour plot is not very different 
from the usual contour. 

• The fraction of interlopers, that is, of galaxies that pass 
the criteria, but are not located in the target redshift range 
[2.7,3.4], is about 12%. Most of these interlopers are fore- 
ground galaxies at z > 2. 

• The efficiency of the selection, defined as the product 
of the completeness by the confirmation rate, is 85%. 

• The contribution of clustering to cosmic variance can 
be estimated for scales below the apparent size of the box at 
redshift ~ 3. For a 100 arcmin 2 field, the relative variance 
is more than twice the Poissonian value. 

• The reconstructed luminosity func tion, provided we ap - 
ply a correction for volume effects as in ISteidel et ail fl995). 



is very similar to the one given in that paper. This shows an 
overall consistency between the mock sample and the actual 
sample. 

Given this overall agreement, we have explored the na- 
ture of LBGs in the model, with special emphasis on three 
critical issues (i) clustering and spatial information, (ii) ex- 
tinction and the optical-to-IR luminosity budget, and (iii) 
merging history and fate in the present universe. We have 
obtained the following results: 

• The amplitude and shape of the projected 2-point cor- 
relation function is in agreement with the data. With our 
choice of the cosmological parameters of the simulation, the 
comparison of the 3D correlation function with the dark 
matter correlation function gives us an estimate of the linear 
bias that slowly decreases on scales from 1 to 20 /i _1 Mpc, 
from b ~ 5 to b ~ 3. So we find a higher bias value than 
previously found in the literature. This is the first estimate 
of the bias for z ~ 3 LBGs from a hybrid model. The LBGs 
are located in haloes with masses ranging from 3 x 1O 11 M0 
to 3 x 10 12 A/ Q , with a median value 1.6 x 1O U M . The po- 
sitions of the LBGs nicely delineate the dense regions of the 
filaments at z ~ 3. 

• The extinction properties of our model LBGs depend 
on the extinction model that has been put in, basically the 
one presente d in the STARDUST mo del of spectrophotomet- 
ry evolution dDevriendt et al]ll999l) . We find a factor ~ 6.2 
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for extinction at the rest-frame w avelength 160 A th at is, 
a value higher than that used bv lSteid el et alj il996t), but 
consis tent with the later analysis of ^A^^iberger^ Steidell 
(2000). The UV to IR luminosity b udget is quite similar 
to the re-construction proposed by lAdelberger fc Steidell 
( 2000) on the basis of the available data. 

• We estimated a set of properties for our mock LBGs. 
The median values are 2.6 x 1O 1O M for stellar masses, 2.6 
kpc for half-light radii, 0.2 Zq for metallicities, and 27 (resp. 
50) M©yr _1 for instantaneous (resp. averaged over 100 Myr) 
SFRs. These median values for stellar masses, radii and 
metallicities are in agreement with the data, although the 
measurement of these quantities from observations is not 
straightforward. Our value for the median SFR seems to be 
too low with respect to observations 3 . We also find that only 
a minority of LBGs (30%) may be starburs ts triggered by 
mergers, contrary to ISomerville et al] <l200l]) . 

• Finally, we are able to assess the future of these z ~ 3 
LBGs. We give examples of their SFR and mass histories 
from this epoch to the present. We are also interested in 
the nature of LBGs as progenitors of current galaxies of 
the Hubble sequence, although morphological classification 
issues are not easy. We recall that in GALICS, our Hubble 
type is given by the B-band bulge-to-disk luminosity ratio. 
The model LBGs observed at z ~ 3 preferentially end up 
in galaxies that are classified as bright, early-type galax- 
ies. More specifically, 77 % of the descendants of our model 
LBGs are ellipticals or lenticulars with Mb < —18, whereas 
the latter population only represents 15 % of the local galax- 
ies with Mb < —18. On the other hand, 35 % of the bright 
ellipticals and lenticulars of the local Universe actually had 
a LBG progenitor at z ~ 3. The other 65 % may have had 
already progenitors at z ~ 3, but with apparent magnitudes 
that are to faint to pass the chosen selection threshold. So 
we cannot say that LBGs "a la Steidel" are the progenitors 
of local ellipticals, but they clearly are progenitors of local 
ellipticals. 

Our model gives a very good basis for making predic- 
tions of observational strategies for forthcoming follow-up 
observations. For instance, we show that only a small frac- 
tion of the LBGs can be detected in a SCUBA follow-up 
or with the SWIRE 4 survey with SIRTF. These objects are 
clearly a key target for ALMA which should be able to de- 
tect all of them at 850 /im or 1.3 mm, at the 0.1 mjy level. 
The strong potential of our model for making follow-up pre- 
dictions through many optical and IR/submm filters can 
be fully exploited by submitting queries to our relational 
database at http://galics.iap.fr 

7.2 Summary of drawbacks 

The main shortcoming of our model for the present study 
comes from the mass resolution of the cosmological simula- 

3 We believe, however, that this underestimate is due both to 
differences in stellar population models used to retrieve the SFR 
and in the definition of the period on which the SFR is com- 
puted. When playing the game of fitting our LBG colours with 
e.g. the PEGASE model iFioc fc Rocca- Volmerangel \l99% . we 
obtain larger values of the SFR, by a factor 2—3. 

4 http : //www . ipac . caltech . edu/ SWIRE/ 



tion we use. Only haloes with 20 or more particles are iden- 
tified by the fof algorithm. This results in a baryonic mass 
limit, below which incompleteness comes in. Galaxies with 
baryonic masses below the threshold are included in identi- 
fied haloes, in which not all the hot gas has cooled, whereas 
galaxies in dark matter structures less massive than our de- 
tection threshold are missing. Fortunately, we have shown 
in GALICS I that incompleteness settles slowly for decreasing 
masses. On its turn, the baryonic mass threshold converts 
into an absolute magnitude limit that we call our formal 
magnitude limit. The estimate of this value at z ~ 3 for 
rest-frame 1600 A absolute magnitude is Mab = —20.3. 
This means that this mass resolution just passes the selec- 
tion criterion put as Mab < —20.2. 

Unfortunately, the mass resolution was shown to "leak" 
on galaxy properties through their merging histories by GAL- 
ICS I. It was not crucial for bright galaxies at z = 0, but it is 
much more significant at z ~ 3. Galaxies that should have 
been missed because they are not massive enough appear 
over the magnitude threshold because they are "immature" . 
Their history is not resolved enough in the past, because of 
the absence of small haloes, and they appear too bright inas- 
much as they are too gassy and form stars too actively. We 
showed that we can recover from this effect by imposing a 
"maturity criterion" : only galaxies with at least a progenitor 
1.1 Gyr ago have to be selected. The other objects have to 
be discarded from our analysis. Such a criterion admittedly 
appears as a fine-tuning to recover the correct number of 
LBGs. It is only at z = 3 that we are 'drowning' in major 
resolution effects, inasmuch as the correction reaches about 
75 %, and we are consequently reaching the limit of this 
specific simulation. A study on a higher-resolution simula- 
tion (Blaizot et al, 2004, in progress) shows that the cho- 
sen maturity criterion becomes transparent at z = 3 in this 
simulation, and that our maturity criterion to by-pass the 
problem in the current simulation is reasonable. 

Another shortcoming appears as the consequence of the 
volume of our cosmological simulation. The comoving side 
of the box 100ft -1 Mpc corresponds to an apparent 1 deg 
at z ~ 3. It does not allow us to address the issue of cosmic 
variance in fields larger than typically a fraction of a deg 2 , 
and clustering on scales larger than a fraction of a degree. 
On scales smaller than that, we have shown in a companion 
paper (momaf) that the technique of cone building from 
this simulation introduces systematic effects that lower the 
correlation function at most by 10 %. 

Finally, we have tried to model the selection criteria as 
accurately as possible, but we may have missed systematic 
effects and subtle observational biases. 



8 CONCLUSIONS 

The LBG selection is a very efficient way to select high- 
redshift galaxies. The z ~ 3 LBGs seem to naturally come 
out of models of hierarchical galaxy formation, as shown by 
the current study. The list of properties that are reproduced 
is quite impressive, and the overall feeling is that we globally 
understand what is going on, both in the selection process 
and in the nature of the target galaxies. However, the actual 
impact of this set of results is somewhat affected by the 
need to correct for the "immaturity" of a large number of 
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z ~ 3 galaxies in the model. This immaturity reflects mass 
resolution in terms of history resolution. For this reason, 
there is no point in attempti ng to go at higher redshift, 
and try to model z ~ 4 LBGs iSteidel et alJll99Sft with the 
outputs of this specific simulation. 

The GALICS model seems to nicely mimic the properties 
of optically-bright, star-forming galaxies at z ~ 3. How- 
ever, this obviously does not mean that other galaxy popula- 
tions in this redshift range that would be selected according 
to other observational criteria, are equally well reproduced. 
More specifically, samples of galaxies may be selected for 
instance on the basis of their red colours, such as the Ex- 
tremely Red Objects. In this prospect, another subtle effect 
of mass resolution may appear in terms of the absence of 
old stellar populations that should have formed in struc- 
tures smaller than the mass threshold. This issue will be 
addressed in forthcoming work. 

After our first two GALICS papers that presented a de- 
tailed exploration of the nature of model galaxies between 
z — and 3, this study is the next milestone in the pre- 
sentation of the full sets of results coming out of the GAL- 
ICS model, as far as it illustrates how mock samples can be 
extracted with elaborated selection criteria. This study is 
complemented by a relational database that can be easily 
queried through a web interface at |htt p : / /galics . iap . f r 
Two forthcoming papers will be more oriented toward pre- 
dictions for deep surveys, in terms of multi-wavelength faint 
galaxy counts (GALICS iv) and 2D correlation functions 
(GALICS V). 

The road map for further developments is rather clear. 
On the one hand, since we are reaching the limit of this spe- 
cific simulation, we have to run simulations with a better 
resolution. Such a work is in progress and seems to nicely 
confirm the results of the current study. On the other hand, 
the logic of mock observing cones should lead us to gener- 
ate mock images with all the instrumental effects, and to 
try to extract mock samples with the same data processing 
pipelines as the actual observations. 
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